Patterns of the Aharonov-Bohm oscillations in graphene nanorings 
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Using extensive tight-binding calculations, we investigate (including the spin) the Aharonov-Bohm 
(AB) effect in monolayer and bilayer trigonal and hexagonal graphene rings with zigzag boundary 
conditions. Unlike the previous literature, we demonstrate the universality of integer (hc/e) and 
half-integer (hc/2e) values for the period of the AB oscillations as a function of the magnetic flux, 
in consonance with the case of mesoscopic metal rings. Odd-even (in the number of Dirac electrons, 
A'^) sawtooth-type patterns relating to the halving of the period have also been found; they are 
more numerous for a monolayer hexagonal ring, compared to the cases of a trigonal and a bilayer 
hexagonal ring. Additional more complicated patterns are also present, depending on the shape 
of the graphene ring. Overall, the AB patterns repeat themselves as a function of N with periods 
proportional to the number of the sides of the rings. 

PACS numbers: 73.23.-b, 73.22.Pr, 73.23.Ra 
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I. INTRODUCTION 

Due to widespread interest in nanoscience and nan- 
otechnology in the past fifteen years, persistent currents 
(PCs) and the Aharonov-Bohm (AB) effect in ring-type 
nanosystems have attracted much attention. Originally, 
PCs and the AB effect were studied theoretically for spin- 
less electrons in the ideal case of strictly one-dimensional 
(zero-width) metallic nanorings threaded by a solenoidal 
magnetic flux.^"^ Subsequently, consideration of spin in 
this ideal case was shown^ to lead to a nontrivial odd- 
even behavior, associated with halving ($o/2 versus $o) 
of the universal AB period and of the corresponding am- 
plitude of the AB oscillations as a function of the applied 
magnetic field B] $o — hc/e is the unit of magnetic flux. 

Recently fabricated carbon-based new materials, like 
carbon nanotubes^ and two-dimensional graphene, pro- 
vide additional opportunities for investigations of PCs 
and the AB effect, with potential future technological 
applications, in ring-type nanodevices. However, in spite 
of the recent extraordinary interest in graphene (start- 
ing with the isolation of a single graphene sheetZ.), only 
a few experimental^'^ and theoretical studies (see, e.g., 
Refs. [l0l - [l3l ) of PCs and the AB effect in graphene nanor- 
ings have appeared in the last couple of years. Surpris- 
ingly these graphene-ring studies have been inconclusive 
regarding the aforementioned odd-even behavior associ- 
ated with the electron spin; at the same time, no regu- 
lar behavior or other pattern of the AB oscillations was 
reported. Moreover, oneii of these publications has con- 
cluded that the odd-even behavior fails to manifest in 
graphene nanorings at all. 

In this letter, based on extensive tight-binding calcu- 
lations, we investigate the AB oscillations for the case 
of trigonal and hexagonal narrow graphene rings termi- 
nating in zigzag edges; for experimental advances in the 
fabrication of graphene samples with well-defined high- 
purity edges, see Ref. 14 Our systematic studies (in 
the size range 1 < < 100 Dirac electrons) reveal clear 
signatures of several well defined patterns (including odd- 
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FIG. 1. (Color online) Diagram of the narrow trigonal and 
hexagonal graphene rings with zigzag boundary conditions 
(for both the inner and outer edges) used in the TB calcu- 
lations, (a) equilateral trigonal ring with a width of three 
rows of carbon atoms (b) hexagonal ring (with edges forming 
concentric regular polygons), with a width of five rows. The 
length unit is the lattice constant a = 0.246 nm. 



even and halved-period behaviors) that can be traced to 
consideration of both the spin degree of freedom and of 
the zigzag boundary conditions obeyed by the graphene 
Dirac electrons. The different conclusion arrived in this 
Letter in comparison with previous publicat ion o ^ ° ' ^ ^ ap- 
pears to be due to the simplified^^ condition (infinite- 
mass boundary condition, which, unlike the zigzag condi- 
tion, cannot describe different crystallographic termina- 
tions and corner geometries in graphene) used in the lat- 
ter, in conjunction with the circular symmetry required 
for obtaining analytic solutions of the continuous Dirac- 
Weyl equation. 



II. PRELIMINARY THEORETICAL 
BACKGROUND 

The spectra of an ideal metallic ring^ (IMR) are very 
regular exhibiting a parabolic dependence on the mag- 
netic fiux $, which is portrayed by the simple analytic 
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where the single-particle angular momentum / takes the 
values I — 0, ±1, ±2, .... This regularity is directly re- 
flected in AB related quantities, such as the presistent 
current I and the total magnetization M, which exhibit 
a periodic behavior as a function of $ with period $o 
(for spinless electrons^) or both $o and $o/2 (when the 
electron spin is considered.—) Indeed one has, 
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where the total energy 
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is given by the sum over all occupied single-particle 
(noninteracting electronsi^) energies; the index a runs 
over spins. The magnetic flux in Eq. ^ is specifled as 
$ — BS, where the area S = ttR^, with R being the 
radius of the ID ideal ring; for advances in the measure- 
ment of small persistent currents and magnetic moments, 
see Ref.|l3(b). 

To determine the single-particle spectrum (the energy 
levels Si{B)) in the tight-binding (TB) calculations for 
the graphene rings, we use the hamiltonian 
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with <> indicating summation over the nearest- neighbor 
sites The hopping matrix element 
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tij = t exp [ j ds ■ A(r) 



where t = 2.7 eV, and rj are the positions of the 
carbon atoms i and j, respectively, and A is the vec- 
tor potential associated with the applied perpendicular 
magnetic fleld B. The diagonalization of the TB hamilto- 
nian [Eq. ^] is implemented with the use of the sparse- 
matrix solver ARPACKi^il In calculating Etot [see Eq. 
©J, only the single-particle TB energies with ei{B) > 
are considered. '^"■^^ We note here that, unlike the contin- 
uous Dirac-Weyl equations, both the K and K' val- 
leys are automatically incorporated in the tight-binding 
treatment of graphene nanorings. 



III. MONOLAYER TRIGONAL RING 

First we analyze TB results for a narrow trigonal 
graphene ring having pure zigzag terminations for both 
the inner and outer edges; see Fig.[T](a). The correspond- 
ing TB spectra are displayed in Fig.^a). Since the con- 
stant magnetic field B is applied across the whole width 




FIG. 2. (Color online) Characteristic TB low-energy spectra 
of the narrow graphene rings with zigzag boundary conditions 
portayed in Fig. [T] (a) trigonal graphene ring, (b) hexago- 
nal graphene ring. The thick black lines indicate the highest 
occupied state for A'' = 15 [in (a)] and A'' = 43 [in (b)] Dirac 
electrons (spin included). Note the three-fold energy bands 
for the trigonal ring in (a) and the six-fold ones for the hexag- 
onal ring in (b). In the case of the trigonal ring (odd number 
of sides), the consecutive three- fold bands are shifted by a 
phase $o/2 with respect of each other; this results to a dou- 
bling of the period of the AB patterns as a function of A'^, 
i.e., a period of twelve instead of six (spin included). In the 
case of the hexagonal-ring spectrum (even number of sides), 
no such shift is present, and the period as a function of A'^ 
remains twelve (spin included). 



of the ring, the magnetic flux is deflned here in an av- 
erage sense, i.e., through the use of an average area S'av 
given by 
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(5'inn + >S'out)/2, 



(6) 



where the indices "inn" and "out" indicate the areas en- 
closed by the inner and outer edges of the ring, respec- 
tively. 

The graphene-ring spectra in Fig. [2ja) arc different 
from the simple spectra in Eq. ([T]), familiar from the 
case of ID metallic rings. ^ Specifically, they are grouped 
in bunches of six levels (see also Ref. TS ), and each 
such bunch contains two three-level units. Naturally, 
this organization will be refiected in the behavior of the 
Aharonov-Bohm oscillations. Indeed, we found that the 
AB oscillations for the magnetization Af(<I>) exhibits an 
overall period of 2x6 = 12 as a function of the electron 
number TV (the factor of 2 resulting from the spin de- 
gree of freedom). Within this period of 12 electrons, we 
find four distinct patterns as a function of $ (see Fig. [3]) , 
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FIG. 3. (Color online) Magnetization as a function of the 
the magnetic flux <& (spin is included). The panels portray 
the four characteristic patterns of the Aharonov-Bohm os- 
cillations associated with the trigonal graphene ring having 
zigzag boundary terminations; see Fig. [TJa). (a) Sawtooth, 
(b) Pinched sawtooth, (c) Asymmetric rounded sawtooth, 
(d) Halved-period sawtooth. 



namely (a) sawtooth, (b) pinched sawtooth, (c) asym- 
metric rounded sawtooth, and (d) halved-period saw- 
tooth. 

The first three patterns [Fig. EJa-c)] exhibit a period 
of $0 as a function of $, while the fourth pattern [Fig. 
^d)] has a halved period $0/2. As aforementioned, the 
halving of the fundamental period $0 was seen earlier in 
studies^ of the AB effect for spinfuU electrons in ideal ID 
metallic rings. In this case, it was described as an odd- 
even effect due to a two-electron alternation as a function 
of N. In contrast, the halving of the fundamental period 
in the case of trigonal graphene nanorings exhibits a six- 
electron period as a function of N, namely for N = 6i + 
No, and only when Nq = 3 (i = 1, 2, . . .). 

Another regular behavior in the AB patterns of trigo- 
nal graphene rings is a constant shift of the ^-dependence 
by ±$0/2 for all electron sizes related by iV = 6i -I- A^Oi 
with A^o = 1,2,..., 6; Nq is kept constant while i runs 
over 1 = 1,2,3,.... For example, the pattern of A^ = 8 is 
the same as that of A^ = 2, but shifted by $o/2, and the 
same holds for the pattern of A^ = 10 relative to that of 
Af = 4, etc. 

Taking consideration of the above, and through inspec- 
tion of magnetization curves in the range 1 < A^ < 100, 
the following summary of the AB patterns can be de- 
duced (i = 1, 2, . . .): 

1. Sawtooth pattern (a) with zero shift: A^ = 12i + 1, 
N = 12i + 2, N = 12i-H 10. 

2. Sawtooth pattern (a) with a $0/2 shift: A^ = 12i -I- 4, 
N ^12i + 7, N = 12i + 8. 

3. Pinched sawtooth pattern (b) with a $0/2 shift: 
A^ = 12i + 5. 

4. Pinched sawtooth pattern (b) with zero shift: A^ = 
12i-|- 11. 




FIG. 4. (Color online) Magnetization as a function of the the 
magnetic flux <& (spin is included). The panels portray six 
patterns of the Aharonov-Bohm oscillations associated with 
the hexagonal graphene ring having zigzag boundary termi- 
nations [see Fig. [TJb)]. The first five patterns [(a)-(e)] corre- 
spond to single-particle states near the middle of the 12-fold 
spectral band (spin included), while the sixth pattern corre- 
sponds to the top state [see Fig. [Hb)]. (a) = 29; shifted 
halved-period sawtooth pattern, (b) N — 30; shifted saw- 
tooth, (c) A'^ — 31; halved-period sawtooth pattern, (d) 
= 32; sawtooth pattern, (e) A'^ = 33; halved-period saw- 
tooth, (f) — 36; rounded sawtooth. The qualitative devel- 
opment in [(a)-(e)] of an odd-even alternation between one- 
period, $0, and halved-period, "l>o/2, sawtooth patterns is 
evident. 



5. Asymmetric rounded sawtooth pattern (c) with zero 
shift: N =12i + 6. 

6. Asymmetric rounded sawtooth pattern (c) with a <&o/2 
shift: N = 12i + 12. 

7. Halved-period sawtooth pattern (d) with zero shift: 
A^ = 12i-|-3. 

8. Halved-period sawtooth pattern (d) with a $0/2 shift: 
N ^12i + 9. 

To summarize, $o/2-oscillations as a function of the 
magnetic flux occur only in cases 7 and 8 above, with 
the latter involving also an overall ^o/2 shift. 



IV. MONOLAYER HEXAGONAL RING 

Next we analyze the AB oscillations in the case of a 
narrow hexagonal graphene ring with zigzag edges [see 
Fig. [ijb)]. The corresponding energy spectrum [see Fig. 
[5Jb)] exhibits again an organization in bands, as was the 
case with the spectra of the trigonal ring. However, each 
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band now contains six, instead of three, single-particle 
levels, and this is clearly connected to the sixfold point- 
group symmetry of the regular hexagon (the three-level 
bands arising also from the threefold symmetry of the 
equilateral triangle). 

Compared to the trigonal-ring spectra, the hexagonal- 
ring spectra are simpler in one way; namely, there is 
no phase shift between two successive sixfold bands [see 
Fig. Elb)], in contrast to the $o/2 shift between succes- 
sive threefold bands for the trigonal rings [see Fig.[5l^a)]. 
The presence (absence) of a $o/2 shift between successive 
bands appears to be a general behavior of the spectra of 
regular-polygon-shaped graphene rings with odd (even) 
number of sides. 

The absence of a shift between consecutive energy 
bands leads to a simplification of the Aharonov-Bohm 
patterns, since it results in a period of 2 x 6 = 12 (avoid- 
ing the doubling to 24) electrons as a function of N. Of 
particulat interest is the fact that, disregarding a po- 
tential shift of ±$o/2, the AB patterns exhibited by the 
magnetization curves (see Fig.^]) display a well developed 
(although apparently not perfect) alternation pattern be- 
tween integer periods ($o) and halved periods ($o/2), 
as long as the highest occupied state lies in the inte- 
rior of the sixfold energy band. The $o/2 period reflects 
the zigzag nature of the interior states (which we term 
W-states to distinguish from the zigzag boundary condi- 
tion) ; examples of W-states are given by the thick black 
lines in Fig. When the Fermi level (highest occupied 
state) coincides with a W-state, $o/2-oscillations occur. 
Note that there are four W-states for the hexagonal ring, 
but only one W-state for the trigonal ring. 

In Fig. 131 we display the magnetization curves for sev- 
eral instances of electrons occupying states in the 12- fold 
band with the number of electrons ranging from A'^ = 25 
to A^ = 36 (the doubling 2x6= 12 is due to consideration 
of the electron spin). In the range 27 < A^ < 34, the mag- 
netization curves exhibit an odd-even effect associated 
with the alternation between a whole-period ($o) saw- 
tooth oscillation and a halved-period ($o/2) sawtooth 
pattern (exhibiting also a halved amplitude); examples of 
this behavior are portrayed in Fig. |4[(a)-(e)]. The two 
cases for AT = 25 and A^ = 26, with the 25th and 26th 
electrons occupying the bottom level of the sixfold band, 
exhibit both a full-period ($o) sawtooth behavior. Fi- 
nally, the two electrons occupying the top level of this 
energy band (corresponding to A^ = 35 and A^ = 36) 
exhibit a dissimilar behavior, with the penultimate one 
{N = 35) having a full-period ($0) sawtooth behavior 
and the ultimate one (A^ = 36) showing a full-period 
($0) rounded-sawtooth behavior [see Fig. SJf)]. Natu- 
rally, the aforementined AB patterns repeat themselves 
with a period of 12 electrons. 

In Fig. [SJ we display illustrative magnetization curves 
for the case of a wider hexagonal ring compared to the 
one in Fig.[ljb) (by a factor of 2.4). From an inspection 
of the patterns in Fig.[5l as well as others not shown here, 
we found that the behavior of the AB oscillations in this 
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FIG. 5. (Color online) Patterns of the Aharonov-Bohm oscil- 
lations associated with a wider [in comparison with Fig. [IJb)] 
hexagonal monolayer graphene ring (with 12 rows of carbon 
atoms at each side) having zigzag boundary terminations, (a) 
= 27; halved-period sawtooth, (b) N = 28; sawtooth. 
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FIG. 6. (Color online) Total energy curves (as a function 
of the magnetic flux $) corresponding to the magnetizations 
in Figs. Ilfb) and UJ^c) [case of the thin monolayer hexago- 
nal graphene ring with zigzag terminations portrayed in Fig. 
mb)]. (a) iV = 30. (b) iV = 31. Observe the doubling of the 
frequency and the halving of the amplitude of the oscillations 
as one goes from = 30 (even) to A = 31 (odd). 



wider ring change only in minor ways. Much wider rings 
are needed to reach a substantial modification in the AB 
behavior. 



V. SIMILARITIES WITH THE IDEAL METAL 
RING 

To gain further insight into the appearance of the odd- 
even AB behavior in graphene nanorings with zigzag ter- 
minations (described in Sees. IIIII and IIV[) . we plot in 
Fig.inithe total energy curves, £'tot($) [see Eq. ([3])], as a 
function of the average magnetic flux $ [see Eq. ([5])] for 
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two characteristic cases; namely for iV = 30 and A'' = 31 
discussed earlier for an hexagonal graphene ring, see Fig. 
Stb) and He). 

A remarkable feature of these total energy curves is the 
almost parabolic (cx dependence on the magnetic flux 
(equivalently the applied magnetic field), which exhibit 
a period $0 for iV = 30 (even) and a half period $o/2 
for TV = 31 (odd). The odd-even sawtooth oscillations of 
the magnetization portrayed in Fig. |4] are a direct conse- 
quence of this parabolic dependence given the definition 
of the magnetization as the derivative of the total energy 
with respect to the magnetic flux [see Eq. ©J. 

We have further examined the total energy curves, 
E™^{^) (not shown here), for the case of an ideal metal- 
lic ring, i.e., using the well known analytic energies of Eq. 
([1]) , and have confirmed that their shape consists of sim- 
ilar parabolic segments exhibiting a $0 or a $o/2 period 
for even or odd N, respectively. 

Naturally, this overall parabolic (oc $^) dependence of 
EI^^{^) could have been anticipated due to the original 
parabolic dependence on $ of the single-particle levels 
ej'^^($) [see Eq. ([ij]. However, for graphene rings with 
zigzag terminations, this result is a surprising one, given 
that the associated single-particle spectrum is much more 
complicated; it further indicates that the corresponding 
graphene single-particle energies [associated with the W- 
states, see Sees. IIIII and lTV] are parabolic on $ to a rather 
large degree. 

We further briefly mention here that in preliminary 
calculations we found that graphene rings with armchair 
edge terminations have, in contrast to those with zigzag 
terminations, single-particle spectra with an almost lin- 
ear dependence on <i>, and thus their AB patterns are 
different (as we will describe in detail elsewhere^^) . 



VI. BILAYER HEXAGONAL RING 

Having addressed the appearance of regular trends in 
the AB oscillations of monolayer graphene nanorings, we 
comment next on possible modifications that arise in as- 
sociated bilayer graphene-ring structures. To this end, 
we consider an hexagonal bilayer ring formed by stacking 
two monolayer rings [resembling the arrangement por- 
trayed in Fig. [TJb)] one on top of the other following the 
Bernal prescription. Due to the Bernal-type coupling be- 
tween the two rings, such narrow bilayer graphene rings 
are analogs of the double-ring configurations considered 
recently in the framework of the Aharonov-Bohm effect 
in mesoscopic metallic devices.^*' 

A charecteristic part of the low-energy TB spectra for 
the bilayer ring is displayed in Fig. [Tja). As was the 
case with the monolayer hexagonal rings, the emergence 
of sixfold energy bands persists also for the case of a 
narrow bilayer hexagonal ring. However, the couplings 
between the layers leads to strong modifications within 
each energy band; namely, the three top energy levels 
are strongly compressed compared to the three bottom 




FIG. 7. (Color online) An example for the case of a narrow 
bilayer hexagonal graphene ring with zigzag terminations in 
all twelve sides. The bilayer ring is built by stacking (Bernal 
stacking) two monolayer hexagonal rings resembling the shape 
in Fig. \V[h). (a) Characteristic part of the spectrum. The 
thick black line denotes the level occupied by the 29th electron 
(spin included). The dashed line (online blue) denotes the 
single W-state here, (b) Corresponding magnetization curve 
(for A'^ = 29 electrons) as a function of the the magnetic flux 
$ (spin is included). 



ones. This results in turn in several more complicated 
profiles for the AB oscillations, an example of which is 
displayed in Fig. ^h). From an inspection of Fig. ^a), 
it is also clear that there is only a single well-formed W- 
state that may serve as a Fermi level [see second level 
from the bottom denoted by a dashed line (online blue)], 
and thus a halved-period sawtooth pattern occurs only 
once within the period of twelve electrons (with the spin 
degeneracy being accounted for). 



VII. CONCLUSIONS 

Using TB calculations and taking into account the 
spin, we have demonstrated the universality of the in- 
teger ($0) and half-integer ($0/2) magnetic-flux periods 
in the Aharonov-Bohm effect in narrow graphene rings 
with zigzag boundary conditions (trigonal and hexago- 
nal shapes were considered in both monolayer and bilayer 
structures). The AB patterns for the monolayer hexago- 
nal rings are dominated by an odd-even (in the electron 
number) alternation of sawtooth-type oscillations with 
$0 and $0/2 periods. Such an odd-even alternation per- 
sists also for trigonal monolayer and hexagonal bilayer 
rings, with a reduced occurrence frequency (related to 
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the number of W-states in each energy band). Addi- 
tional patterns of higher complexity are also prominent, 
depending on the structure of the graphene ring. All AB 
patterns repeat themselves as a function of N with peri- 
ods relating to the point-group symmetry of the geomet- 
rical shape of the rings. Our findings, which contrast 
with the results of recent literature on the subject (see, 
e.g., Refs. [13 and [III), provide the impetus for experi- 



mental probing of the AB effects in the graphene systems 
explored in this paper. 
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